Possibilities of inversion of satellite third-order gravitational tensor onto gravity anomalies: a case study for central Europe

Wed 04 Sep 2019 09:55:16 AEST

We investigate a numerical performance of four different schemes applied to a regional recovery of the gravity anomalies from the third-order gravitational tensor components (assumed to be observable in the future) synthetized at the satellite altitude of 200 km above the mean sphere. The first approach is based on applying a regional inversion without modelling the far-zone contribution or long-wavelength support. In the second approach we separate integral formulas into two parts, that is, the effects of the third-order disturbing tensor data within near and far zones. Whereas the far-zone contribution is evaluated by using existing global geopotential model (GGM) with spectral weights given by truncation error coefficients, the near-zone contribution is solved by applying a regional inversion. We then extend this approach for a smoothing procedure, in which we remove the gravitational contributions of the topographic-isostatic and atmospheric masses. Finally, we apply the remove-compute-restore (r-c-r) scheme in order to reduce the far-zone contribution by subtracting the reference (long-wavelength) gravity field, which is computed for maximum degree 80. We apply these four numerical schemes to a regional recovery of the gravity anomalies from individual components of the third-order gravitational tensor as well as from their combinations, while applying two different levels of a white noise. We validated our results with respect to gravity anomalies evaluated at the mean sphere from EGM2008 up to the degree 250. Not surprisingly, better fit in terms of standard deviation (STD) was attained using lower level of noise. The worst results were gained applying classical approach, STD values of our solution from T_zzz are 1.705 mGal (noise value with a standard deviation 0.01 × 10⁻¹⁵m⁻¹s⁻²) and 2.005 mGal (noise value with a standard deviation 0.05 × 10⁻¹⁵m⁻¹s⁻²), while the superior from r-c-r up to the degree 80, STD fit of gravity anomalies from T_zzz with respect to the same counterpart from EGM2008 is 0.510 mGal (noise value with a standard deviation 0.01 × 10⁻¹⁵m⁻¹s⁻²) and 1.190 mGal (noise value with a standard deviation 0.05 × 10⁻¹⁵m⁻¹s⁻²).

Comparison of spherical cap and rectangular harmonic analysis of airborne vector gravity data for high-resolution (1.5 km) local geopotential field models over Tanzania

Thu 13 Apr 2023 09:38:45 AEST

This study examines local geopotential field modeling over a mountainous region in Tanzania using vector airborne gravity data. We use the adjusted spherical cap and rectangular harmonic analyses. Both methods are based on expansion of gravitational potential into a series of orthogonal harmonic basis functions of local support in such a way that the expansion coefficients are determined by gravity observations. All three components of gravity vector are simultaneously inverted to derive the geopotential coefficients. In order to evaluate the accuracy of the local models, independent checkpoints are selected within the study region and around its boundary and the computed gravity vectors are compared with the independent gravity observations. The results show an excellent agreement with root mean square error (RMSE) of < 1.6 mGal over the study area. On the contrary, the RMSEs of global geopotential models against the checkpoints data are 7 mGal for the models up to the maximum degree of 2190 (a resolution of ∼9.1 km) and 5 mGal to 5399 (∼3.7 km). Our local models are significantly more accurate than the state-of-the-art global models and fully exploit the airborne vector data with the measurement error of ∼1 mGal. We also present the regional models constrained only by radial (vertical) or by lateral (horizontal) gravity observations. Those models are considerably less accurate than the one from 3D gravity data inversion. Lastly, the regional models are validated against topography data. It is found that the gravity-topography correlation is 0.8–0.9 at 100 km, 0.5 at 20 km, and higher than the correlations of the global models at all frequencies. The gravity-topography admittances estimated from our regional models indicate ∼130 mGal/km and imply the effective density of 2500 kg/m3 for topographical mass.